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Abstract 

New solutions are found for the non-relativistic hydrodynamical equa- 
tions. These solutions describe expanding matter with a Gaussian density 
profile. In the simplest case, thermal equilibrium is maintained without 
any interaction, the energy is conserved, and the process is isentropic. 
More general solutions are also obtained that describe explosions driven 
by heat production, or contraction of the matter caused by energy loss. 

Introduction. The equations of hydrodynamics correspond to local conser- 
vation of some charges as well as energy and momentum. The equations are 
scale-invariant, hence can be applied to phenomenological description of physi- 
cal phenomena from collisions of heavy nuclei to collisions of galaxies. Recently, 
a lot of experimental and theoretical effort went into the exploration of hydro- 
dynamical behaviour of strongly interacting hadronic matter in non-relativistic 
as well as in relativistic heavy ion collisions, see for example refs. Due 
to the non-linear nature of the equations of hydrodynamics, exact solutions of 
these equations are rarely found. In ref. Q an exact solution of hydrodynamics 
of expanding fireballs was found 20 years ago. The purpose of this Letter is 
to present and analyze a new, exact solution of the non-relativistic hydrody- 
namical equations, with a generalization to heat production or loss (e.g. due to 
radiation). We hope that the results presented herewith may be utilized to ac- 
cess analytically the time-evolution of the hydrodynamically behaving strongly 
interacting matter as probed by non-relativistic heavy ion collisions |,|]. The 
results are, however, rather general in nature and they can be applied to any 
physical phenomena where the non-relativistic hydrodynamical description is 
valid. 

Adiabatic expansion. Consider a hydrodynamical system described by the 
continuity equation, the Euler equation and the local energy conservation: 

^ + V(vn) = 0, (1) 
d \ VP 



^+vV V = , (2) 
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|£+V(ve) = -PVv. (3) 
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We close this set of equations with the help of an ideal gas equation of state, 

e{r,t) = ^P(r,t), (4) 
P(r, t) = nir,t)T{t). (5) 

These equations are solved by 

TV f \ 

m . R^t) . RMt), v(r,t) = ||r, (7) 

= + + = (8) 

where mq = Rq/t is the characteristic expansion rate. This is a new solution of 
the non-relativistic hydrodynamical equations ||l|. The flow field has the same 
distribution as in the solutions found by Zimanyi, Bondorf and Garpman 20 
years ago however the temperature and the density profiles are different. 

Now consider the phase-space distribution function / belonging to this hy- 
drodynamical solution. This distribution has a locally Maxwell-Boltzmann 
(MB) shape everywhere, all the time: 

ff n ( (p-mv(r, t)f \ N 
/(r, p, t) = Cexp ^ —— , C = 



2R^{t) 2mT{t) y (47r2mToi?^)3/2 • 

(9) 

Suppose that the distribution is created at i = to by an instantaneous source. 
It solves the following coUisionless Boltzmann equation for t > 0: 

I + W) / ^ Sir, p, t) . Cexp (-^ - 1^^) Sit ^ to). (10) 

Since the system is considered after to, the normalization coefficient C is con- 
stant because Tit)R^it) — T^Rq = const. Solutions to the coUisionless Boltz- 
mann equation are not uncommon, see e.g. however the cited solutions 
generate deviations from the MB shape of the local momentum distributions, 
hence they are not solutions of the hydrodynamical equations. 

In contrast, assuming here a special initial condition, spherical symmetry and 
a radial initial flow pattern, i.e. eq. (|To|), we find a solution of the coUisionless 
Boltzmann equation that maintains a MB form everywhere for all times. In 
this sense, our solution is really a solution of hydrodynamics. Physically this 
corresponds to a gas that first had been thermalized in equilibrium. Later, 
the gas was frozen out in an idealized break-up dA, t — to- Our local phase- 
space distribution function /(r, p, t) maintains its locally thermalized shape 
without any collisions, due to the special choice of the initial conditions and the 
invariance of the Gaussian shape under convolution. 
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We find a non-vanishing pressure maintained without colhsions; the inter- 
pretation of this result is that any wall or bubble inserted into this expanding 
Knudsen gas would feel a pressure that arises due to the random, locally dis- 
ordered motion of the free-streaming particles in any part of space. Such a 
pressure arising from a coUisionless gas of photons is well-known in cosmology 
as a source of gravity, see ref. ||^. 

The entropy density can be evaluated from kinetic theory, using (|^): 

By introducing the "effective volume" Vq = (27r)'^/'^i?g and integrating the 
entropy density over space, we find that the total entropy is 

S{t) = 5ideai + ^iV = const. (12) 

As it was shown in jlO| , one can modify the thermodynamical definition of the 
entropy by adding terms linear in extensives, without changing the thermody- 
namics. Above, we got the interesting result that the entropy is almost the 
same as S'idoai, the entropy of an ideal gas at temperature Tq in volume Vg 
— the difference is the extensive quantity |iV, hence the thermodynamics of 
the system considered is the same as that of an ideal gas. Note that the total 
entropy is independent of time. It is easy to verify from eqs. (|^, |l^) that the 
process is locally isentropic: 

|^ + V(v.) 0. (13) 

Without rescattering and other final state interactions after particle emis- 
sion, the momentum spectrum is independent of time: 

^r(p) = mip, t) = ^J^^ exp (-^) , ^ To + mul (M) 

It is worthwhile to evaluate the total energy in local disordered motion (heat 
energy, denoted by ii'hcat), the total energy in ordered motion (flow energy, 
denoted by inflow) and the total energy i?tot- One obtains that 

3 3 
^^tot = -NT^ = const, E^hoat = -NT{t), Efio„ = £^tot - £^hcat, (15) 

where the time dependence of the local temperature is given by eqs. (0-||). 

The two-particle spectra A^2(Pi, P2, ^i, ^2) are also independent of time if 
ti,t2 > tQ. The two particle distribution function can be written in terms of the 
emission function S in the following way: 

tl t2 

^2(pi, P2, h, t2) = J dt\J rfVi J dt'2 J d\2 (5(ri, pi, t[)S{r2, P2, Q 

— 00 —00 

± ^(n, K/2, t\)S{v2, K/2, t'2) cos [q(ri - V2)] , 
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where K = (pi + P2)/2, q = Pi — P2 and the "+" sign refers to bosons, "-"to 
fermions, if the final state Coulomb and other interactions are negligible. S 
vanishes for t[ 2 > so the the time integration is independent of the upper 
bound ti,2 as long as ti^2 > to- Therefore N2 is independent of ti and t2- 

Because both the single and the two-particle spectra are independent of time, 
the two-particle correlation function must also be independent of time: 

I \ ^2(Pl, P2) , -i?2q2 

C2 Pi, P2 = -iTTi — Vl^TT — T = 1 + e *^ ■ (16) 

A^l(Pl)A^2(P2) 

So the radius parameter of the correlation function, i?* is time independent. 
We can calculate it at toi the time of particle production. Let uq = Rq/t. The 
result 

1111 



i?J i?^ Rq \ To ) ^ ^ 

is a generalization of the result in ref. Q for arbitrary values of uq. 

More general solutions. The above presented solution of the non-relativistic 
hydrodynamical equations can be generalized in many ways. For instance, a 
straightforward way would be to introduce an inhomogeneous temperature pro- 
file lO] . Such a temperature profile is present in a known analytic solution, see 
ref. |p[. Here we will investigate the possibility in our model to describe the 
effects of some local heat production or heat loss. A simple model is introduced 
to mimic the effects of heat production by chemical or nuclear reactions, or 
the cooling of the system by radiation that decreases the local energy density. 
The microscopic details of these processes are not sought for and the sources 
of heat production, or the radiated energies are not part of the system under 
consideration, so heat production of radiation changes the 'total' energy. 

Let us introduce an additional term into the energy balance equation (^). 

de 3 

- + V(ve) + FVv = -j(i)n(r, t)Tit). (18) 

This new term is the simplest possible model of heat production {j{t) > 0) or 
energy loss e.g. due to radiation {j{t) < 0). The heat loss or heat production is 
assumed to be proportional to the local internal energy. 

Now we have — in place of eqs. — a more general Gaussian solution of 
the continuity equation: 

m = To^, R^t) = v(r,i) = ^r, (19) 

where h is so far an unknown function. This parameterization satisfies the 
boundary conditions for the temperature and the radius if 1^9(^0) = ^(^o) = 1- 

By substituting the expressions of the energy density (^, the pressure (^ 
and the particle density (^ into the modified energy equation (^, and using 
the parameterization (^^, we get the following differential equation for h: 

hit) = hit)jit). (20) 
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Note that adding source terms to the energy balance equation results is a de- 
viation from the isentropic expansion or contraction (]l3|). From (pi]), (^, ( p^ 
and (pO|), the local entropy production is 

|^ + V(v.) = ^jn. (21) 



The solution of eq. (pOl) is 



/.(i) = e^^o''''^''\ (22) 



The function (^(t) can be determined from the Euler equation, 

0{tMt)~-^'{t) = —^h{t). (23) 
Supposing j{t) = const ^ 0, equation (|23) has two exponential solutions: 



^(t) = e^^--^*-*"), h(t) = eH*-*o)^ ^- = ± /l^. (24) 
In terms of temperature, radius and velocity, the solution reads as: 



T{t) ^ Toe5J(*-*«), R\t) = i?2e^j(*-to)^ ^(j., t) = ^r. (25) 

In case of j > 0, we have an exponential expansion and warming up driven 
by an external energy source. The opposite case is the exponential contraction 
of the matter, for j < 0. It may be caused by the continuous emission of energy 
by radiation. The flow field is independent of time in both cases. 

Algorithm to generate new solutions. One can generate even infinitely many 
new analytical solutions to the equations of non-relativistic hydrodynamics with 
energy producing or radiative processes, using the following method: 

1) Assume a functional form for (p{t) and choose {TQ,m, Rq) so that 

0{to)--^Hto) = ^{to) = 1. (26) 



2) Determine the function h(t) from eq. (p3[). 

3) Find the energy source function j{t) from eq. (pO|). 

This way, a solution of eqs. (|lH|, |l^) is generated. The only requirement for 
consistency arises from T{t) > 0, which results in h{t)/Lp{t) > 0. The solution 
is given by eq. ( p^ ) with ip(t), h{t) and j{t) generated by steps l)-3). 

Summary. We found a new class of solutions for the non-relativistic hydro- 
dynamical equations, describing a spherically expanding Knudsen gas with time 
dependent but location independent temperature. 

Then we incorporated into our formalism the possibility to describe the 
effects of the emission or absorption of heat. In case of the simplest analytic 
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solution, the system expands or contracts exponentially, while the flow field 
becomes independent of time in this special case. 

Finally, we presented an algorithm that generates infinitely many new ana- 
lytical solutions of the hydrodynamical equations with energy sources. 
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